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Abstract 

The complete fermionic two-loop contributions to the prediction for the W-boson 
mass from muon decay in the electroweak Standard Model are evaluated exactly, i.e. 
no expansion in the top-quark and the Higgs-boson mass is made. The result for the 
W-boson mass is compared with the previous result of an expansion up to next-to- 
leading order in the top-quark mass. The predictions are found to agree with each 
other within about 5 MeV. A simple parametrization of the new result is presented, 
approximating the full result to better than 0.4 MeV for Mh < 1 TeV. 
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The prediction of the W-boson mass, Mw, in terms of the Z-boson mass, Mz, the Fermi 
constant, G^, and the fine structure constant, a, is one of the most important quantities for 
testing the electroweak Standard Model (SM) and its extensions with high precision. This 
relation is derived from muon decay, as the Fermi constant is defined in terms of the muon 
lifetime, r^, according to 




1 + Ag), (1) 



with F{x) = 1 — 8x — 12x^lna; + 8x^ — x'^. By convention, the QED corrections within 
the Fermi Model, Aq, are included in this defining equation for G^. The one- loop result for 
Aq pp, which has already been known for several decades, has recently been supplemented by 
the two-loop correction [2] . The tree- level W propagator effects giving rise to the (numerically 
insignificant) term 3m^/ (5M^) in eq. (P) are conventionally also included in the definition 
of Gfj,, although they do not belong to the Fermi Model prediction. 

Comparing the prediction for the muon lifetime within the SM with eq. yields the 
relation 

where the radiative corrections are summarized in the quantity Ar This relation can be 
used for deriving the prediction of Myv within the SM or extensions of it, to be confronted 
with the experimental result for Mw- At present, the W-boson mass is measured with an 
accuracy of 5 x 10~^, = 80.419 ± 0.038 GeV jlj. The experimental precision on Mw 

will be further improved with the data taken at LEP2 in its final year of running, and at the 
upgraded Tevatron and the LHC, where an error of (5Mw = 15 MeV can be expected 
At a high-luminosity linear collider running in a low-energy mode at the W"'"W~ threshold, 
a reduction of the experimental error down to 5Mw = 6 MeV can be envisaged This 
offers the prospect for highly sensitive tests of the electroweak theory provided that the 
accuracy of the theoretical prediction matches the experimental precision. 

The one-loop result for Ar within the SM [3] can be decomposed as (with = 1 — 

Ar(") = Aa- ^Ap + Ar^^Mu), (3) 

where the leading fermion-loop contributions Aa and Ap, arising from the charge and 
mixing- angle renormalization, are separated out, while the remainder part Arj-em contains 
in particular the dependence on the Higgs-boson mass, Mr. The QED-induced shift in the 
fine structure constant, Aa, contains large logarithms of light-fermion masses. The leading 
contribution to the p parameter from the top/bottom weak isospin doublet, Ap, gives rise 
to a term with a quadratic dependence on the top-quark mass, rrit [H]. 

Beyond the one-loop order, resummations of the leading one-loop contributions Aa and 
Ap are known [9J. They correctly take into account the terms of the form (Aa)^, (Ap)^, 
(AaAp), and (AaArrom) at the two-loop level and the leading powers in Aa to all orders. 
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While the QCD corrections to Ar are known at 0{aas) ^U] and 0{aa^) jTI], only partial 
results are available up to now for the electroweak two-loop contributions. They have been 
obtained using expansions for asymptotically large values of rut [Tallin] and Mr [H]. The 
terms derived by expanding in the top-quark mass of O^C^mf) [12] and 0(G^mtM|) [T^ 
were found to be numerically sizeable. The 0{G'jjnlM'^) term, involving three different mass 
scales, has been obtained by two separate expansions in the regions Mw, Mz, Mh -C rrit and 
My^^Mi -C mt,MH and by an interpolation between the two expansions. This formally 
next-to-leading order term turned out to be of a magnitude similar to that of the formally 
leading term of 0{G'j^mf), entering with the same sign. Its inclusion (both for Mw and the 
effective mixing angle sin^ ^eff) had important consequences on the indirect constraints on 
the Higgs-boson mass derived from the SM fit to the precision data. 

Consequently, a more complete calculation of electroweak two-loop effects appears de- 
sirable, where no expansion in rrit or Mh is made. As a first step in this direction, exact 
results have been obtained for the Higgs-mass dependence (e.g. the quantity Mw,subtr(^H) = 
Mw(Mh) — Mw(Mh = 65 GeV)) of the fermionic two-loop corrections to the precision ob- 
servables fT^. They have been compared with the results of expanding up to (9((j^m^M|) 
specifically analysing the effects of the rrit expansion, and good agreement has been found 

Beyond the two-loop order, complete results for the pure fermion-loop corrections (i.e. 
contributions containing n fermion loops at n-loop order) have recently been obtained up 
to four- loop order [Q. These results contain in particular the contributions of the leading 
powers in Aa as well as the ones in Ap and the mixed terms. 

In the present paper, all fermionic two-loop corrections to Ar are calculated exactly, i.e. 
without an expansion in the top-quark or the Higgs-boson mass. These are all two-loop dia- 
grams contributing to the muon decay amplitude and containing at least one closed fermion 
loop (except the pure QED corrections already contained in the Fermi model result, see 
eq. (^). Figure H displays some typical examples. The considered class of diagrams includes 
the potentially large corrections both from the top/bottom doublet and from contributions 
proportional to iYjf and N^^, where A^if is the number of light fermions (a partial result for 
the light-fermion contributions has been obtained in Ref. JH])- The results presented here 
improve on the previous results of an expansion in rrit up to next-to-leading order in 
containing the full dependence on rrit as well as the complete light-fermion contributions at 
the two-loop order, while in Ref. higher-order corrections from light fermions have only 
been taken into account via a resummation of the one-loop light-fermion contribution. 




Figure 1: Examples for types of fermionic two- loop diagrams contributing to muon decay. 
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In the following, we briefly outline the main features of the calculation. After extracting 
the IR-divergent QED corrections that are already contained in the Fermi model QED factor 
(a detailed description of how this is done will be given in a forthcoming paper), the generic 
diagrams contributing to muon decay can be reduced to vacuum-type diagrams, since the 
masses of the external particles and the momentum transfer are negligible. The on-shell 
renormalization of the gauge-boson masses, on the other hand, requires the evaluation of 
two-loop two-point functions with non-zero external momentum, which is more involved 
from a technical point of view regarding the tensor structure and the evaluation of the 
scalar integrals. It should be noted that this complication cannot be avoided by performing 
the calculation within another renormalization scheme (the MS scheme, for instance), since 
ultimately one is interested in the relation between the physical parameters Mw, Mi-, a, G^, 
rather than between their MS counterparts. For this reason we have decided to use the on- 
shell renormalization scheme everywhere in our calculation, i.e. we use physical parameters 
throughout (alternatively one could of course do the calculation in a different renormalization 
scheme, with formal parameters, and perform the transition to the physical parameters in a 
second step) . If not otherwise stated, we use the conventions of Ref . [1^] . 

In our calculation we have made use of some computer-algebra tools. The package 
FeynArts [201 ™^ applied to generate the Feynman amplitudes and counterterm contri- 
butions. The program TwoCalc [21] was applied for the algebraic evaluation of these am- 
plitudes, which were reduced, by means of two-loop tensor-integral decompositions, to a set 
of standard scalar integrals. The calculation was carried out in a general gauge, which 
allowed us to test the gauge-parameter independence at the algebraic level as a highly non- 
trivial check. For the evaluation of the scalar one-loop integrals and the two-loop vacuum 
integrals we have used analytical results as given in Ref. [221 ? while the two-loop two-point 
integrals with non- vanishing external momentum have been evaluated numerically using one- 
dimensional integral representations with elementary functions . These allow a very fast 
calculation of the integrals for general mass configurations. 

Since we are using Dimensional Regularization |24,.;25j in our calculation, a careful treat- 
ment of the Dirac algebra in D dimensions involving 75 is necessary. While a naively anti- 
commuting 75 can safely be applied for all two-loop two-point contributions (for a discussion, 
see e.g. the first paper of Ref. [121) most of the two-loop vertex- and box- type diagrams, 
this is not the case for the two-loop vertex diagrams containing a triangle subgraph, shown 
in Fig. 121 For these graphs, a naively anticommuting 75, although respecting the Ward 
identities, would lead to an incorrect result. This is due to an inconsistent evaluation of 
the trace of 75 together with four Dirac matrices, which in four dimensions is given by 
Tr {757'^7''7''7'^} = Aie^^^" , while applying the naively anticommuting 75 in D dimensions 
would yield zero for this trace. In order to calculate this type of diagrams, we have first 
evaluated the triangle subgraph with the mathematically consistent definition of 75 in D 
dimensions according to Refs. [231201 (here we made use of the package Tracer [27] for 
checking). After adding appropriate counterterms, which are necessary to restore the Ward 
identities, the result differs from the result obtained using a naively anticommuting 75 only 
in terms proportional to the totally antisymmetric tensor e^^p'^ , Inserting the latter contri- 
bution into the two-loop diagrams, we find that the second loop gives a finite contribution. 
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so that it can be evaluated in four dimensions without further comphcations.^ The fermion 
hne appearing in the second loop also yields an e-tensor contribution, which results, after 
contraction with the e-tensor from the triangle subgraph, in a non- vanishing contribution to 
the result for Ar. 




Figure 2: Two- loop vertex diagrams containing a triangle subgraph, which require a careful treat- 
ment of 75 in D dimensions. 

As mentioned above, we perform the renormalization within the on-shell scheme. It 
involves a one-loop subrenormalization of the Faddeev-Popov ghost sector of the theory, 
which is associated with the gauge-fixing part. The gauge- fixing part is kept invariant under 
renormalization. For technical convenience, we manage this by a renormalization of the gauge 
parameters in such a way that it precisely cancels the renormalization of the parameters and 
fields in the gauge-fixing Lagrangian.^ To this end we have allowed two different bare gauge 
parameters for both W and Z, ^J^'^ and and also mixing gauge parameters, ^"^^ and 

. The renormalized parameters comply with the gauge, with one free gauge parameter 
for each gauge boson. With this prescription no counterterm contributions arise from the 
gauge-fixing sector. Starting at the two-loop level, counterterm contributions from the ghost 
sector have to be taken into account in the calculation of physical amplitudes. They follow 
from the variation of the gauge-fixing terms under infinitesimal gauge transformations. 
We have derived all the counterterms arising from the ghost sector (extending the results 
of Ref. [201 to a general i?^ gauge) and implemented them into the program FeynArts. In 
this way we could verify the finiteness of individual (gauge-parameter-dependent) building 
blocks (e.g. the W- and the Z-boson self-energy) as a further check of the calculation. 

Concerning the mass renormalization of unstable particles, from two-loop order on it 
makes a difference whether the mass is defined according to the real part of the complex 
pole of the S matrix, 

= M^ - iMT, (4) 

or according to the pole of the real part of the propagator. In eq. (0} A4 denotes the complex 
pole of the S matrix and M, F the corresponding mass and width of the unstable particle. 
We use the symbol M for the real pole. 

-'^For recent discussions of practical ways of treating 75 in higher-order calculations, see also Refs. |29II28| . 

^An alternative way of achieving that the gauge-fixing sector does not give rise to counterterm contribu- 
tions would have been to add the gauge-fixing part to the Lagrangian only after renormalization, in which 
case the renormalized gauge transformations would have to be used. 
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In the context of the present calculation, these considerations are relevant to the re- 
normalization of the gauge-boson masses, M-w and Mz- The two- loop mass counterterms 
according to the definition of the mass as the real part of the complex pole are given by 

«w,(2) = Re{Si'p,(AC)}-M4_,i,iZ«+Im{E?;;,,(A4)}lm{E?;,i,(A4)}, (5) 

z 

+ Im {S|5;)(M|)} Im {e|5,)(M|)} , (6) 

where St,(i), St,(2) denote the transverse parts of the one-loop and two-loop self-energies 
(the terms from subloop renormalization are understood to be contained in the two-loop 
self-energies), and means the derivative of the one-loop self-energy with respect to the 

external momentum squared. Field renormalization constants are indicated as 5Z^ . The 
relations to the mass counterterms according to the real-pole definition, 5M^ 5M| (2), 

are given by 

5M^,(2) = 5M^,(2) + Im{s^|^)(M^)}lm{sWi)(M^)}, (7) 

5M|,(2) = 5M|,(2)+Im{s|^'i)(M|)}lm{s^^,)(M|)}. (8) 

It can easily be checked by direct computation that the terms in eqs. (|Z|), (|H1) by which the two 
definitions differ are gauge-parameter-dependent. Thus it is obvious that at least one of the 
two prescriptions leads to a gauge- dependent mass definition. While the problem of a proper 
definition of unstable particles in gauge theories has already been addressed many times in 
the literature jSI], it should be noted that in the present calculation two- loop contributions of 
the type leading to a non-zero (and gauge-parameter-dependent) difference between the two 
kinds of mass renormalization methods are for the first time fully included in a computation 
of a physical observable in the Standard Model. Explicitly, these are contributions from 
light fermions and bosonic loops evaluated in a general gauge. In the previous results for 
Mw, incorporating terms up to O^G'j^m^M^) [T^j and Mn-dependent fermionic terms [T3] . 
the contribution Im |sfp (-^)(M^)| Im |St,(i)(M^)| was zero, making thus a strict distinction 
between the two mass definitions unnecessary at the considered order. 

Since our result has been obtained within a general gauge, we can investigate the 
issue of whether the mass renormalization is gauge-parameter-independent by explicit com- 
putation. In particular, the two-loop counterterm to the weak mixing angle, 5sw,(2), ought 
to be gauge-parameter-independent since sw is a physical observable (note, however, that 
the same argument does not hold for the mass counterterms of eq. © and eq. (jH]); see e.g. 
Ref. 122] for a discussion). We find that 5sw,(2) is only gauge-parameter-independent with 
the definition of the gauge-boson masses according to the complex pole, while the real-pole 
definition for the masses leads to a gauge-parameter-dependent result for 5sw,(2)- This result 
is in accordance with what is expected from S-matrix theory, in which the complex pole is 
a gauge- invariant quantity |ST| . 

We have thus adopted the complex-pole definition as given in eq. (0) and eq. (0). Us- 
ing this mass definition leads to a Breit-Wigner parametrization of the resonance line 
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shape with a constant decay width. Experimentally the gauge-boson masses are deter- 
mined using a Breit-Wigner function with a running (energy-dependent) width. Con- 
necting the latter prescription with the theoretical prediction involves the approximation 
Im |s^^^^(s)| ^ srw,z/^w,z, which is vahd for the fermionic contributions to the W- and Z- 
boson self-energies at one-loop order. As usual, Fw.z denote the W- and Z-boson widths. As 
a consequence of the different Breit-Wigner parametrizations, there is a numerical difference 
between the mass parameters corresponding to the definition used in the experimental deter- 
mination (denoted as Mw, Mi henceforth) and the mass parameters in our calculation, Mw, 
Mz. The shift between these parameters is given by ^w,z = ^w,z + r^.z/(2Afw,z)- 
Since M-w and Mz enter on a different footing in our computation — Mz is an experimental 
input parameter, while Mw is calculated — in order to evaluate the mass shifts we use the 
experimental value for the Z-boson width, Fz = 2.944 ± 0.0024 GeV and the theoretical 
value for the W-boson width, which is given by F-w = 3G^M^/(2^/27r)(l + 2as/(37r)) in 
sufficiently good approximation. This results in Mz ~ Mz + 34.1 MeV and in the mass 
shifts Mw ~ Mw + 27.4 MeV and Mw ~ Mw + 27.0 MeV for Mw = 80.4 GeV and 
Mw = 80.2 GeV, respectively. =^ 

We now turn to the numerical discussion of our result for Ar. It should be noted that 
our definition of Ar according to eq. is based on the expanded form (1 + Ar) with 
Ar = Ar*^") + Ar*^" ^ + . . . rather than on the resummed form 1/(1 — Ar), indicating a 
resummation of leading one-loop contributions. The terms consistently taken into account at 
two-loop order with such a resummation are explicitly contained in our two-loop contribution 
to Ar. The result for Ar contains the following contributions 

Ar = Ar^") + Ar^""^) + Ar^""') + Ar^^'"') + Ar^^'"'^ (9) 

where Ar*^"-* is the one-loop result, eq. 0, Ar^""*"-* and Ar*^""=^ are the two-loop ^Hj and 
three-loop ^T] QCD corrections, while Ar'-^*" ^ is the new electroweak two-loop result. The 
notation [NiO^) symbolizes the contribution of all diagrams containing one fermion loop, 
where iVf stands both for the top/bottom contribution and for all light-fermion species. The 
term Ar^^^f contains the pure fermion-loop contributions in two-loop order. Since the pure 
fermion-loop contributions in three- and four-loop order have been found to be numerically 
small, as a consequence of accidental numerical cancellations, with a net effect of only about 
1 MeV in Mw (using the real-pole definition of the gauge-boson masses) ^7j, we have not 
included them here. 

In Fig. Elthe different contributions to Ar are shown as a function of Mh. Here Mw is kept 
fixed at its experimental central value, Mw = 80.419 GeV, and mt = 174.3 GeV is used. 
The effects of the QCD corrections, of the two-loop corrections induced by a resummation 

■^The difference in Fw according to the way it is calculated, through the tree-level result parametrized with 
a, or the improved Born result parametrized with G^, or the improved Born result including QCD corrections 
(which is the one we used), is formally of higher order (i.e. beyond Oia^)) in the calculation of Afyv- Its 
numerical effect is nevertheless not completely negligible; it changes the shift in Mw by about —2.9 MeV if 
the tree- level result for Fw parametrized with a is used and by about —1.4 MeV if the parametrization 
of the Born width (without QCD corrections) is employed. 
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of Aa, and of the purely electroweak fermionic two-loop corrections are shown separately. 
The purely electroweak two-loop contributions are sizeable and amount to about 10% of 
the one-loop result. We have compared the Higgs-mass dependence of Ar with the result 
previously obtained in Ref. J3] and found perfect agreement. 
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Figure 3: Different contributions to Ar as a function of Mh- The one-loop contribution, Ar'"\ is 
supplemented by the two-loop and three-loop QCD corrections, Arq^j-, = Ar'""=^ -|- Ar'""s), and 
the fermionic electroweak two-loop contributions, Ar*^"^) = Ar''^'"^^ -|- Ar'^f^"^^. For comparison, 
the effect of the two-loop corrections induced by a resummation of Act, Arj^^\ is shown separately. 

The prediction for Mw is obtained from the input parameters by solving eq. (j21). Since Ar 
itself depends on Mw this is technically done using an iterative procedure. The prediction 
for Mw based on the results of eq. is shown in Fig. 0] as a function of Mh for = 
174.3 ± 5.1 GeV jS3 and Aa = 0.05954 ± 0.00065 The current experimental value, 

= 80.419 ± 0.038 GeV [H, and the experimental 95% C.L. lower bound on Mh (Mh = 
107.9 GeV 36 ) from the direct search are also indicated. The plot shows the well-known 
preference for a light Higgs boson within the SM. Confronting the theoretical prediction 
(allowing a variation of rrit, which at present dominates the theoretical uncertainty, and Aa 
within la) with the la region of and the 95% C.L. lower bound on Mh, only a rather 
small region in the plot (corresponding to 107.9 GeV < Mh ^ 140 GeV) matches all three 
constraints. 

We have compared our results with those of an expansion for asymptotically large values 
of rrit up to 0(G^m^M|) [iHllHZl- The results are shown in Table [T] for different values of 
Mh. For the input parameters the values of Ref. [THj have been chosen, i.e. rrit = 175 GeV, 
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Figure 4: The SM prediction for Mw as a function of Mh for mt = 174.3 ±5.1 GeV is compared 
with the current experimental value, M-^^ = 80.419 ± 0.038 GeV [1], and the experimental 95% 
C.L. lower bound on the Higgs-boson mass, Mh = 107.9 GeV [HS] . 



Mz = 91.1863 GeV, Aa = 0.0594, a^{Mz) = 0.118. Relatively good agreement is found, 
with maximal deviations in M-w of about 5 MeV. If we had chosen a different parametrization 
of Fw in the above calculation of the shift between the masses corresponding to the fixed and 
the running width definition, a somewhat larger deviation to the result of Ref. ^13j would 
have been obtained. 

The deviations in the last column of Table ^ can of course not be attributed exclusively 
to differences in the two-loop top-quark and light-fermion contributions, because the results 
also differ by a slightly different treatment of those higher-order terms that are not yet under 
control, such as purely bosonic two-loop contributions and effects from scheme dependence. 
A detailed discussion of those differences and of the remaining theoretical uncertainties from 
unknown higher-order corrections will be given in a forthcoming publication. 

Following Ref. we also provide a simple numerical parametrization of our result for 
Mw- It is given by 

Mw = - ci dH - C5 dH^ + cq dH^ - ca da + Cg dt - cj dH dt - C4 da,, (10) 

where 

= ^HiooGiv j ' = lirlaGiv j = 0:05924-^' = Tn^-^' ^''^ 

and Mz = 91.1875 GeV 4J has been used. For the coefficients ci, . . . , C7 we have obtained via 
a least squares fit = 80.3755 GeV, ci = 0.05613, ca = 1.081, C3 = 0.5235, C4 = 0.0763, 
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Mh/ GeV 


Mw/ GeV 


M^PV GeV 


AMw/ MeV 


65 


80.3985 


80.4039 


-5.4 


100 


80.3759 


80.3805 


-4.6 


300 


80.3039 


80.3061 


-2.2 


600 


80.2509 


80.2521 


-1.2 


1000 


80.2122 


80.2129 


-0.7 



Table 1: The two- loop result for Mw based on eq. Q is compared with the results of an 
expansion in rrit up to O^Gf^m'^M^) OEZl, M^^'^. The last column indicates the difference 
between the two results. 

C5 = 0.00936, cq = 0.000546, C7 = 0.00573. The parametrization of eq. (fTUj) approximates 
our full resuh for Mw within 0.4 MeV for 65 GeV < Mh < 1 TeV. 

In summary, we have evaluated the complete fermionic two-loop contributions to the 
W-boson mass within the electroweak Standard Model. Our result improves on previous 
results as it does not involve any approximations in the top-quark and the Higgs-boson mass 
and also contains the contributions of all light fermions in the Standard Model. Within 
our calculation we have defined the gauge-boson masses according to the complex pole of 
the S matrix, which ensures the gauge-parameter independence of the mass definition. We 
have provided a simple numerical parametrization of our result, which approximates the full 
result with sufficient accuracy for all values of Mh up to 1 TeV. In comparison with the 
previous result obtained for Mw by an expansion for asymptotically large values in rrit up 
to next-to- leading order we find slightly lower values of Mw, sharpening thus the tendency 
towards a light Higgs boson within the Standard Model. 
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